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Bubble nucleation has been studied on lattices using phenomenological Langevin equa- 
tions. Recently there have been theoretical motivations for using these equations. These 
studies also conclude that the simple Langevin description requires some modification. 
We study bubble nucleation on a lattice and determine effects of the modified Langevin 
equations. 

1. Introduction 

The effective potential determines ttie vacuum state of a quantum field theory.^ 
When radiative corrections are included or, if one constructs the finite tempera- 
ture version of field theory the (effective)potential may have two minimum, one of 
which is lower than the other. The 'higher' minimum is referred to as the false or 
metastable vacuum. The global minimum is the true vacuum state of the theory. 
The existence of these metastable vacua have interesting consequences in the early 
Universe. If at early times the field settles in the false vacuum it may 'decay' to the 
true vacuum either by quantum tunneling or thermal hopping. The transition region 
is spherically symmetric in coordinate space. Bubbles of the stable phase appear in 
the metastable phase. In relativistic field theory a formalsim for calculating decay 
rates from the unstable to the stable vacuum has been developed.^ Since potentials 
which are non-linear in the fields are inherent in these problems, it is usually not 
always possible to find analytic solutions to the decay rate equations. Numerical 
methods are employed in order to avoid too many analytic approximations. For 
example, in one detailed numerical study, a scalar field in 1 + 1 dimensions in con- 
tact with a thermal bath was modeled using a phenomenological Langevin equation 
with a field independent, white noise driving term.^ There have been attempts to 
derive a Langevin equation containing fluctuation and dissipation terms from purely 
field theory considerations.'' In these works the authors considered different models 
for the thermal bath, but the common result obtained was the possiblity that the 
thermal noise in the Langevin equation could be colored and depend on the field. In 
this work we investigate the effects of these non-linear fiuctuation and dissipation 
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terms numerically. 

2. Numerical Method 

We study nucleation by evolving the equations of motion of the field on a lattice 
using a staggerd leap-frog algorithm. We work in 1 + 1 dimensions and assume that 
the decay is thermally driven. The field is coupled to the thermal bath through a 
stochastic noise term, ^{x,t). 

- - VFW., 0) ^ = -VW + t)G(0(x, t)) (1) 

where F{(j){x,t)) — G{(f){x,t)) — 1 for additive noise and, F{(f>{x,t)) — (j){x,t)'^ and 
G{(j){x,t)) = (f)(x,t) for multiplicative noise. The general shape of the potential is 
show in Fig. 0(a). The noise is assumed to be Gaussian, in terms of the distribution 
of sizes of the fluctuations, and white, in that the fluctuations are uncorrelated in 
space and time. We assume the noise and viscosity terms are related through the 
fluctuation-dissipation theorem.^ 

< ^{x, t)^{x',t') >= 2Tri5{x ~ x')8{t - t'). (2) 

In later work a generalized fluctuation-dissipation theorem for the multiplicative 
noise will be discussed. We start with the field in the false vacuum, V{(t)) — 0, and 
evolve the equation of motion on the lattice until the field acquires its true vacuum 
value at enough neighboring space points to avoid counting as bubbles fluctuations 
that re-collapse to the false vacuum value. 




(a) (b) 

Fig. 1. (a)The potential V{<j>) = <j>^ — a<j>^ +\4>^ vs (j>: (b)The correlation < (f>{x)(f>{x') > vs \x — x'\ 
for initial conditions at different temperatures 

3. Results 

We report on some observations and preliminary results of our study. A large 
number of nucleation times has to be obtained in order to make a good estimate of 
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the nuclcation rate. Furthermore, the nucleation times obtained from the simulation 
have to be suitably interpreted. We found that if one starts in the false vacuum 
with a uniform field configuration or with one which is random and uncorrelated 
then there is a long waiting time before the any bubble is nucleated. This delay 
is interpreted as the time taken for the coupled field to acquire the short range 
correlation that is needed for a bubble to nucleate. Nucleation of bubbles occurs 
only when the field at some point in space reaches the true vacuum and pulls along 
neighbouring points. A quenching technique, where the equations are evolved in 
an unbroken potential until the short distance correlation, Fig.0(b), is obtained, 
can be used to eliminate the delay time. Fig. ^(a) and (b) show the nucleation 
times recorded for 5000 bubbles using separately the additive and multiplicative 
noise from (|l|). We used a — 0.74, A — 0.1, and 77 = 1. At the time of quenching 
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Fig. 2. Nucleation times at (a) T=1.6 and (b) T=l.l. Ther sharper peaks are the additive noise 
cases. 

we reset the time to zero. We see that the distributions are not symmetric. This 
is expected for the distribution of waiting times for a decay process. We also see 
that the for the multiplicative noise case the nucleation times are larger and more 
broadly distributed. Although the size of the thermal fluctuations is multiplied 
by the field value, the nonlinear dissipation term has an overall larger effect and 
retards the decay in comparison to the additive noise case. If one obtains non- 
linear Langevin type equations when studying realistic models like the Electroweak 
effective potential, for example, the implication for the time scale of nucleation will 
require a similar careful study. Other aspects of these non-linear models currently 
under study include the effects on bubble wall expansion and the spatial distribution 
of critical fluctuations. 
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